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O '. Abstract 

o 

We evaluate N dependences of correlation functions in the bosonic part of the IIB 
matrix model by the Monte Carlo method. We also evaluate those in two sorts of 
regularized Schild models and find that the N dependences are different from those in 
the matrix model. In particular, the distribution of the eigenvalues are logarithmically 
divergent in the regularized Schild model when g'^N is fixed. 
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§ '. 1 Introduction 

O ■ The IIB matrix model is the zero volume limit [@] of a ten-dimensional super Yang-Mills 
■ theory, which is expected to be a non-perturbative formulation of superstring theory. The 
^ . gauge fields in the IIB matrix model are represented by x Hermitian matrices. The 
D ! distribution of eigenvalues of those matrices is interpreted as our space-time.Q Then, how 
^ ' the IIB matrix model settles the space-time structure, e.g., space-time dimension, is one of 

the most important problems. 
^ ■ We need to search for the true vacuum of the IIB matrix model to solve the above 

problem. When A^ = 2, we can actually integrate the fermions and investigate analytically 
the model with the bosonic degrees of freedom . However, it is almost impossible to follow 
the same procedure when N > 2. So we should take other ways to analyze the IIB matrix 
model. It is natural to think that the IIB matrix model in the large- A^ limit is described 
by a field theory because it has infinite numbers of degrees of freedom. We expect that we 
may be able to determine the true vacuum of the IIB matrix model with the techniques of 
field theories if we could find such a field theory. We have paid attention to the fact that 
there is a one-to-one correspondence between the action of an A^ x A^ matrix model and that 
of the field theory on a non-commutative periodic lattice (NCPL) with N x N sites, which 
we denote as T^^ i, |, 0, i. Since the commutation relations of coordinates on NCPL are 
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^This model represents an open space-time since the eigenvalues of the matrices can take values in the 
infinite region. Some studies of a closed space-time with periodic boundary conditions were done with the 
unitary matrices instead H. 
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proportional to one may think that the matrix model in the large- iV limit has a one- 
to-one correspondence to a field theory on the continuous torus. Since the action of a naive 
large- limit of the IIB matrix model is the supersymmetric Schild action on a torus 0, it 
is believed that the finite- IIB matrix model is a regularized supersymmetric Schild model 
|T0| . For this reason, the supersymmetric Schild model on a torus is a plausible candidate 
for the large- iV limit of the IIB matrix model. 

However, we pointed out a possibility that the large- behavior of the IIB matrix model 



and that of the supersymmetric Schild model with the momentum cut-off are different |11 



And the differences between the bosonic matrix model and the lattice-regularized Schild 



model were discussed in [Q. In our previous paper |TT[], we proposed a method to estimate 
the large- A^ dependences of correlation functions of the IIB matrix model kinematically. We 
also applied this method to the supersymmetric Schild model with the momentum cut-off 
and found that it has a different large- A^ dependence from the IIB matrix model because its 
action has a different overall factor of A^. 

In this paper, we confirm our conjecture by calculating correlation functions using Monte 
Carlo simulation. As for the bosonic matrix model, our Monte Carlo results are, of course, 
consistent with those in Ref . . We also calculate some correlation functions in two kinds 
of regularized Schild model, one of which uses the lattice regularization and the other in- 
troduces a naive momentum cut-off. We find that the bosonic matrix model and both of 
the regularized Schild model have different large- A^ dependences. Comparison between the 
large- A^ dependences of the bosonic matrix model and those of the lattice regularized Schild 
model were studied before in Ref.[|l2] and our results are consistent with theirs. This paper 
is organized as follows. In the next section we show the numerical results by the Monte 
Carlo simulation of the bosonic matrix model and in section 3, we consider the regularized 
Schild models and present the results by the Monte Carlo simulations. The final section is 
devoted to summary and discussion. 



2 The bosonic matrix model 



In Ref. [jrT| , we proposed a method to estimate the large- A^ dependences of correlation func- 
tions of the IIB matrix model kinematically. Here we restrict our discussions on the bosonic 
part of the IIB matrix model. The action of the bosonic matrix model is given by 

^M = -T\Tr[A^At, (1) 

where are N x N Hermitian matrices and /U, z/ = 1, ■ ■ ■ , 10.0 In this case our kinematical 
study leads to the result that any correlation functions behave as 

lTrF(A^'))<0(ArO), (2) 

when A^ goes to infinity keeping g^N fixed. This is the same statement as 'tHooft argued 
perturbatively 0] and it was shown perturbatively in Ref.[0, while we do not use any 



perturbative arguments but only kinematical naive counting. We should note that the naive 
^Partition function was evaluated to be finite 
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Figure 1: The N dependences of the correlation functions. In the left figure, A represents 
(pTr [A^, A'^]'^) / (10(A^^ — 1)), which can be analytically evaluated to be 1. Other marks, ■, 

□, •, represent {j^TT{A'^)^)/{g^N), {j^ Ti A"^) /{gN^), {j^Ti (A^A-'f) /{g^N), respectively. 
In the right figure, the marks, ▲, ■, T, •, *, □, represent {j^Ti (A'^Y) / {5g^N2)^ 
Tr {A^A^)^)/{5g^Nl), (i Tr (A^r) /{20g^N^), Tr {A^A'^An^)/{5g^Nl), 

{j.^:tA\A^'A^'y) /(S^^A^I), Tr {A^'A''APA^APA'')) /{5g^N^), respectively. 



counting does not lead to eq. (0) with the usual matrix forms. For example, when we calculate 
the dependence of ((1/A^) Tr (A^)^) naively, we find 

^Tr(A^)A = /l {^A^liK))l^O{N'), (3) 

' \ i,j,k,l=l I 

since the summation is over terms. In general, we have 
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^TrF(A^)^ <0(iV^-i), (4) 

where F{x) is a polynomial of degree p. The situation is changed, however, when we map 
the IIB matrix model to the field theory on NCPL P, ^, ^ We can rewrite any matrix 
models to field theories on NCPL because the commutation relations of the plain waves on 
NCPL coincide with U(N) algebra |^. Correlation functions are rewritten by 

^TrF(A'^)^ = <^^5^F.(A'^(a))|), (5) 

where has the same functional form as F{-) but the products of the arguments are 

replaced by the star products on NCPL, a are N"^ points on NCPL, and hence eq.(^) is at 
most 0{N^) because the summation over a is the summation over iV^ points. 

We calculate some correlation functions using Monte Carlo methods and the results are 
shown in Fig.|I[ The results in the left figure of Fig.^ correspond to the previous calculations 



in Ref.p|. Since {{I / g^)TI[A^' , A'^f) / {IQ{N^ - 1)) is evaluated analytically to be 1, it was 
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calculated to check on the validity of our calculations and it gives perfectly the expected value 
of 1. The right figure of Fig.|l| shows the correlation functions of all possible independent 
combinations of homogeneous polynomials of degree six with rotational invariance and a 
particular one of degree eight, ((1/A^) Tr(yl^)^). We calculate these correlation functions to 
show more evidences for eq.(^. When we fit those data to a + hN'^ by the least square 
method, we find that c is negative and c = —2 fits well [0. Fig.|l] shows that all of the 



correlation functions are convergent in the large- limit, which strongly supports eq.(g). 

3 The regularized Schild models 

It is natural to think that the large- limit of the IIB matrix model is the supersymmetric 
Schild model 10|. This will be easy to understand when we rewrite the bosonic matrix 
model to the field theoretical form on NCPL as follows, 

= Tr [A^ A'^f = ^ E [^'(-)' ^^(^)]* ' (6) 

where 

[A{a),B{a% = A{a) ^ B {a) - B {a) ^ A{a) , 
A{a)i.B{a) = A{a) e-'^ ^ ^ B{a). 

In the large- iV limit, if we neglect naively the non-leading terms of eq.(|^) becomes 

and the last equation is exactly the Schild action on the torus. The space where the fields 
live in is the two-dimensional torus because the momenta of the plane waves in finite N are 
discrete |TI1. Thus the supersymmetric Schild model is a possible candidate for the action 



of the IIB matrix model in the large- iV limit. So we calculate correlation functions using 
some regularized Schild model by Monte Carlo simulation to see whether they agree with 
the results of the bosonic matrix model or not. We adopt two kinds of regularizations for 
the Schild model, one is the lattice regularization and the other is the momentum cut-off. 
The field theoretical description of the finite- matrix model is a field theory on NCPL with 
N X N sites, and the commutation relation of the coordinates is proportional to 1/N. So 
it is natural to think that the finite but large- matrix model approximately agrees with 
the Schild model regularized by the commutative lattice with N x N sites. On the other 
hand, the leading term of the field theoretical form of the finite- A^ matrix model action in 
momentum representation exactly agrees with the Schild model action with the momentum 
cut-off. Then we adopt these two regularizations to compare with the finite- A^ matrix model. 
The actions of these regularized Schild model are as follows in the momentum representation: 

N ^° 

= 4^ E E '^(^ - ^' - k, k) (8) 

fi,u=l m,n,k 
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where m, n, k are (mi,m2), (ni,n2), {ki,k2), respectively and run from — (A^ — l)/2 to 
(N - l)/2, and 

5cut-ofr(m, n) = — m x n, 

oiattice(m: n) = — € sm sm , 

latticed , ; V ^ / V ^ / 

where an extra factor has been multiphed in order to compare with the matrix model 
clearly |11|. Since the fields live in the commutative space here, A^{a)A'^{a) = A''{a)A^{a), 
we calculate {{l/N'^)^{A'^{a))'')/{g^N)^/'^, k = 1,2,3,4, to compare with the correlations 
functions of the bosonic matrix model in Fig.|l]. Also we calculate {{l/g'^N)Y^{dA^^ x 
dA")'^) I (10(A^^ ~ 1)) in order to check on the validity of our calculations, since it is evaluated 
analytically to be 1. The results of the correlation functions of these two regularized models 
are given in Fig.|^. The correlation functions do not seem to be convergent so far, or it looks 
divergent in the large- A^ limit. In particular, when we fit the data of the distribution of the 
eigenvalues to a + hN'^ by the least square method, we find c ~ 0. Hence we put 

^alog(A^) + 6, (a, 6 : constants) (9) 

and find that this fits well. All these results are completely different from those of the bosonic 
matrix model, and we can conclude that the bosonic matrix model is not the regularization 




of the Schild model. These results and conclusion are consistent with an early work |]T2 
where they compared the bosonic matrix model and the lattice Schild model. 



4 Summary and Discussion 

We have calculated correlation functions of the bosonic matrix model and two sorts of reg- 
ularized Schild model by Monte Carlo methods to investigate the large- A^ limits of those 
models. We have the following two conclusions: 

(i) The Monte Carlo results of the bosonic matrix model strongly supports eq.(^, which 
says that any correlation functions such as ((1/A^) TrF(A)) is at most 0(A^°) when A^ goes 
to infinity with g^N fixed to some constant. 

(ii) The A^ dependences of correlation functions of the bosonic matrix model do not agree 
with neither of those of regularized Schild models. 

Let us consider the first statement. We see that the upper bound of the A^ dependences 
of the correlation functions can be calculated by naive counting. In fact, there are some ways 
to count A^ dependence naively and the strongest constraint is eq.(|D. Let us demonstrate 
it with ((1/A^) Tr (A^)^) as an example. We have at least three different ways to express it 
by components of matrices as follows: 



(11) 
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Figure 2: The dependences of correlation functions. The left figure is for the cut- 
off regularization, while the right one is for the lattice regularization. In both cases, A 
represents {-^Y^{dA^{(y) x dA^{o)y) / {l{i{N'^ — 1)), which can be analytically evaluated 

to be 1. The marks, *, ■, ., □ represent (7^ E (1^ E(^'(^))')/(^?'^^), 

Ylii.^^ / {9"^ ^) ^ respectively. In fact each expectation value 

is divided by some constant numerical factor, respectively, so that all results can be seen in 
single frame at once. 



(^m-n^n^-m-k^k) ) ) (12) 
\m,n,k / 

where m = (mi,m2),n = (ni,n2),k = (/ci,/c2). The r.h.s. of eq.(|TUp is a usual matrix 
representation, eq.(|ri|) is the coordinate one on NCPL and eq.(|12|) is its momentum one. 
The N dependences by naive counting are at most 0{N^), 0{N^), and 0{N^), respectively. 
In general, the upper bound of dependences of ((1/A^) Tt A^) by naive counting are given 
by 0(A^^), 0{N'^), and ©(A^^*-^"^-*) from the matrix representation, the NCPL coordinate 
one, and its momentum one, respectively. The NCPL coordinate description, eq. (pA]) , gives 
the strongest constraint on the upper bound of large- dependence and they agree with 
eq.(0). We did not take a perturbative approach |T^, [1^ but only the kinematical one, and 



we expect that our arguments are applicable to the supersymmetric one [p|] . 

Next let us discuss the second statement. We should notice that the distributions of the 
eigenvalues, {g'^N)~^^'^{{l/N'^)'^(A{a))'^), diverge logarithmically in the regularized Schild 
models, while it will be a constant in the large- A^ limit of the bosonic matrix model. In the 
IIB matrix model, ((1/A^^) ^(^(c))^), which is expected to be 0{g^N), is interpreted as 



the extent of the universe, while the fixed g'^N will be proportional to a'^ |17], so that the 
origin of the huge ratio should be explained. On the other hand, if the extent of the universe 
diverges as in the regularized Schild model, or we may fix g'^N^~'^ ~ a'^ {k > 0) fl^ in the 
bosonic matrix model, there would be no problem of the huge ratio. From this point of view, 
it is important to consider the possibility that whether one can build the string theory from 



the supersymmetric Schild model instead of the IIB matrix model as in Ref. [|r^ 
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